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We study the interplay of polaronic effect and superconductivity in transport through molecular 
Josephson junctions. The tunneling rates of electrons are dominated by vibronic replicas of the 
superconducting gap, which show up as prominent features in the differential conductance for the 
dc and ac current. For relatively large molecule-lead coupling, a features that appears when the 
Josephson frequency matches the vibron frequency can be identified with an over-the-gap structure 
observed by Marchenkov et al. [Nat. Nanotech. 2, 481 (2007)]. However, we are more concerned 
with the weak-coupling limit, where resonant tunneling through the molecular level dominates. We 
find that certain features involving both Andreev reflection and vibron emission show an unusual 
shift of the bias voltage V at their maximum with the gate voltage V g as V ~ (2/3) V g . Moreover, 
due to the polaronic effect, the ac Josephson current shows a phase shift of -k when the bias eV 
is increased by one vibronic energy quantum hu v . This distinctive even-odd effect is explained in 
terms of the different sign of the coupling to vibrons of electrons and of Andreev-reflected holes. 

PACS numbers: 74.78.Na, 74.25. fc, 74.45. +c, 74.50.+r 



I. INTRODUCTION 

Electron transport through quantum dots embedded 
in Josephson junctions is attracting increasing interest!^! 
Rich phenomena^HSl arise due to the competition of super- 
conductivity and strong interactions. A crucial feature of 
molecular Josephson junctions (MJJs) is the interplay of 
the ac Josephson effect and molecular vibrations. Signa- 
tures in the transport can be expected when the Joseph- 
son frequency matches a vibration frequency. Since the 
Josephson frequency can be precisely controlled by the 
bias voltage, such signatures could form the basis of a 
precise molecular spectroscopy. This idea has recently 
been explored by Marchenkov et aZ.J^who have measured 
the differential conductance of a Nb dimer in a super- 
conducting break junction. An over-the-gap structure 
consisting of a series of peaks was ascribed to the excita- 
tion of vibrational modes by the ac Josephson current,^ 
but no microscopic description was provided. Theoretical 
treatments of MJJs have mainly focused on the dc cur- 
rent either for weak molecule-lead coupli ng a t vanishing 
biaa^ or for weak electron- vibron couplingSIHl so that the 
microscopic understanding of spectroscopic signatures in 
transport through MJJs is still far from complete. How- 
ever, such an understanding is the prerequisite for de- 
veloping the dc and ac Josephson currents in MJJs into 
spectroscopic tools. 

The purpose of this paper is to investigate the trans- 
port properties of a biased MJJ within a microscopic 
model. We focus on MJJs with strong electron-vibron 
coupling. In this limit, electrons dress with a vibron 
cloud, forming polarons. A minimal model of a single- 



orbital molecule between two s-wavc superconducting 
leads is considered. Electrons in the molecular orbital 
are coupled to a vibrational mode of frequency u) v . We 
will see that signatures in transport do not only occur 
when the Josephson frequency ujj equals the vibration 
frequency u> v but also when hwj, hu) v , and the super- 
conducting gap in the leads satisfy certain other simple 
rational relations. 

The remainder of this paper is organized as follows: In 
Sec. [H] we present our model and discuss the theoretical 
approach. In Sec. |III| we then present and discuss our 
results for the dc and ac Josephson effects, followed by a 
summary in Sec. |IV| 



II. MODEL AND METHOD 

Our model Hamiltonian reads 

H= 2_j H a + H nlo \ + HT, (1) 

a=L,R 

where the first term describes the left (L) and right (R) 
BCS superconducting leads, 

CakCT cSj ko .C Q kcr +^ (A* Cak-fC^-k^ +H.C.) , (2) 
kcr kcr 

with superconducting order parameter A Q . c a \^ a (c^ ko .) 
annihilates (creates) an electron of wave vector k and 
spin a in lead a. 
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The molecule with vibration degree of freedom is rep- 
resented by 



H mo \ — 



' £ m d\d a + ftwyCL^a + \(aJ + a) } d\d a , (3) 



where e m is the molecular energy level, d a (d\) is the 
annihilation (creation) operator of a spin-cr electron in 
the molecular orbital, a (a^) is the vibron annihilation 
(creation) operator, and A is the electron-vibron cou- 
pling strength. We neglect the Coulomb interaction in 
the molecule, which is justified if the charging energy is 
small compared to the coupling to the leadsPSl The role 
of the Coulomb interaction in MJJs has recently been 
reviewed in Ref. ^ 

The tunneling between the molecule and the leads is 
described by 



Hi 



aker 

+ H.c, 



2eV a 



(4) 



where <f> a is the initial phase of the superconducting or- 
der parameter at time t = 0, V a is the voltage in lead a, 
and t a d is the tunneling matrix element. In the follow- 
ing, we choose (j) a = 0, Vl = 0, and Vr = — V, where V 
is the voltage drop across the junction. For symmetric 
capacitances between the molecule and the leads, as as- 
sumed here, the molecular energy level is then given by 
£ m = £o ~ eV/2. 

To go beyond perturbative approaches] 3 * 5 * 8 ! we employ 
the unitary Lang-Firsov transformatiorP^HII] to diagonal- 



ize Hn 



The transformed Hamiltonian reads 



H = H n 



E 



so 



eV 



d\d a + huj v a) a + Ht, (5) 



where the molecular energy level is shifted to e = 
Sq — X 2 /huj v by the polaron binding energy. To simplify 
notation, we now take e$ to denote the shifted level. In 
principle, the electron-electron interaction is also renor- 
malized by the transformation but we neglect this shift 
together with the bare Coulomb interaction. We do not 
expect the on-site interaction to qualitatively change our 
results, which are concerned with transport outside of the 
Coulomb-blockade regime. Ht has the same form as Ht, 
except that the tunneling matrix elements are dressed by 
the polaronic effect as t a d = tudX, where 



X 



exp 



A 



{a) - a) 



(6) 



is the polaron-shift operator. 

The transport properties are o btaine d by the non- 
equilibrium-Green-function method! 12 * 14 * In the Nambu 
representation, we introduce the contour-ordered Green 
function G{t,t') = -i{T c ip{t)^ (t')) with the contour- 
ordering directive T c and ip — (df,d^) T . The particle 



current through the lead a yHHIsl 

I a (t) = ^Re| dt^v {<j z [G<{tM)K(h,t) 

+ G r (t,i 1 )£<(i 1 ,i)]}, (7) 

where the trace is over Nambu space, <r 2 is a Pauli ma- 
trix, and S° and are, respectively, the advanced and 
lesser self-energies due to the coupling to lead a. The ad- 
vanced self-energy is related to the retarded one through 

The current measured in, say, the left lead consists of 
the particle current /t(t) plus the displacement current 
due to the formation of image charges. This contribution 
vanishes in the stationary state but must be taken into 
account in the time-dependent case to ensure gauge in- 
variance and current conservation. This requires proper 
partitioning of the displacement current.^ Since we as- 
sume symmetric coupling, the displacement currents are 
symmetric, in which case the measured c urrent equals 
the symmetrized current I = (II — 7r)/2.E2EU Due to 
the ac Josephson effect, the current can be expanded as 
I(t) = J2 n Ine mu}jt , where the Josephson frequency is 
uij = 2eV/h and /„ = (Il,u ~ ^R,n)/2. The Fourier 
components I a , n are 

J m 

+ G-„ jm (e) S< m (e)] }, (8) 
where G m<n (e) = G„_ m (e + muj) and 



G n {e) 



dt- dt' e tt{t - t] e- mu,jt G(t,t') (9) 



oo T Ja 



and analogously for the self-energies. The current can 
be decomposed in to dis sipative (I®) an d nondissipative 
(1%) contributions' 15 ^ 

I{t) = Io + ^2 {I® cos nwjt + In sin nuijt) , (10) 

n>0 

where 1% = Re (I_„ + I n ) and 1% = Im (J_„ - I n ). 

Before turning to transport properties, we show that 
the molecule-lead coupling is drastically modified by the 
polaronic effect. We employ the standard decoupling ap- 
proximation, which assumes that averages of products of 
two polaron-shift operators can be taken out of el ectron ic 
Green functions and evaluated in equilibrium! 12 * 19 ^^ 
This approximation is valid if hu> v or A is large compared 
to the molecule-lead coupling. It is known that beyond its 
range of quantitative reliability, the approximation pre- 
dicts resonant features that are too broad and too low, 
but that are centered at the same energy as obtained 
from a more advanced treatment.^ In contrast to previ- 
ous studies for normal leads, we encounter not only the 
correlation function (X(t-t')X^) but also (X(t-t')X). 
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FIG. 1: Imaginary part of the first diagonal element of the 
retarded self-energy, — Im(E£)n, which represents scatter- 
ing of electrons between the molecule and lead L, for var- 
ious electron-vibron coupling strengths A. We have taken 
hw v = 8 A and Y — 0.2hu) v . With increasing A, multiple 
vibronic replicas of the superconducting gap edges appear. 



These functions are given b} |12 | 13 | 20 J2 ^ 



i e 



(11) 



(X(t)X) = (X^t)X^)=J2(-^Lie- m ^ t , (12) 



where 



Li = e -(A/^) 2 (2iV„+i) 



exp 



2k B T 



h(v), (13) 



with the modified Besscl function of the argument 
r\ = 2(X/Huj v ) 2 yj N V (N V + 1). N v is the average vi- 
bron number at temperature T determined by the Bosc- 
Einstein distribution function. The correlation functions 
can be evaluated an alytically. While the normal one has 
the well-known forrr l 12 ' 13 ' 20 ' 2 ! ! 

(X(i)Xt) = exp^p-e-^C^ + l) 

+ ^(1-6-"*)]), (14) 



with g = X/huj v , the anomalous one is 
(X(t)X) 



exp(-.g 2 [(l + e-^^ + l) 
+ ^(l + e^*)]). (15) 



The two correlation functions differ by a phase shift of 
half a vibration period. This phase shift will turn out to 
be crucial for the polaron effect on the Josephson current. 
Its origin can be traced back to the factor of (—1)' under 
the sum in Eq. (12 1. This factor stems from Andreev re- 
flection: An electron tunneling out of the molecule trans- 
mutes into a hole before it tunnels back in, which couples 
to the vibron with the opposite sign, as seen by inspecting 
the last term in Eq. 

Under the decoupling approximation, we obtain the 
retarded and lesser self-energies due to the coupling to 
superconducting lead a as 

OO s 



2 \ i-'a-mnltj V c i 
1 

2i 



l - (^[p< m jii](er) - m^mnhM))} (is) 



and 



E(-d 

I— — OO 



] iJ (e + lw,,),(17) 



respectively, where = e ± Z/icau, and i,j = 1,2 are 
Nambu indices. Here, S Q is the self-energy in the absence 
of electron-vibron coupling. Taking the wide-band limit 
and making use of our choice o f vani shing initial phases 
<f> a at time t = 0, S Q is given bjff^ 1 



^L:m»( E ) 



1 -- 



. (18) 



^.;mn( e ) 



u in // e m+l/2) 
^m,n+l/?i?(& ' A 



-1/2J 



£m-l/2 



E L;mn( e ) = * F £ <*mn/(e m ) p L (e m ) Aj . 



1 



*m,n-l /?fl(e m +l/2) g m f* 2 



S mn /(em+1/2) Pfl(em+l/2) #m,n-l /(W1/2) PR^m+1/2) i m f* 



(19) 
(20) 
(21) 



S m ,n+\ /(em-l/2) Pr{^w.-1/2) 1 8mn f(e m -l/2) PR(e m -l/2) I' 

m — 1/2 / 

I 

where e m = e + mhuij and /(e) = (e e /' CsT + 1) _1 is the Fermi-Dirac distribution function, where the cher 
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ical potential in the left lead is independent of the bias 
voltage due to our assumption of Vl = and has been 
absorbed into e. The shifts in the arguments e m ±i/2 of 



the Fermi functions in Eq. (21 1 account for the potential 
difference between the leads. T a = 27r\t a d\ 2 Pa describes 
the coupling of the molecule to lead a, where p^ is the 
density of states in lead a in the normal state. In the 
wide-band limit, r Q is independent of energy. Also, for 
symmetric coupling we have Tl = = T. The expres- 
sions for P a and p a read 



/3a (e) = 



for A Q > |e|, 



for A Q < |e|, 



p a (e) = 9(\e\-A a ) 



respectively. Furthermore, 



H[F](u) =\pj de 



F(e) 



(22) 
(23) 

(24) 



is the Hilbert transform of the function F, where P 
denotes the principal value. The above results are 
valid for fast vibron relaxation so that the averages of 
bosonic operators can be taken in equilibrium. Fig- 
ure [T] shows — Im(S£)n for various electron-vibron cou- 
pling strengths A. Without electron-vibron coupling, 
— Im(S£)n exhibits the superconducting gap of the 
left lead. In the presence of electron-vibron coupling, 
— Im(S£)ii develops vibronic replicas of the gap edges, 
separated by integer multiples of huj v , which open in- 
elastic tr anspo rt channels beyond the usual Andreev 
reflection! 15 ! 16 ! They enable electrons with energies above 
the superconducting gap to undergo Andreev reflection 
under emission or absorption of vibrons. 
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FIG. 2: (Color online) Density plot of the dc differential con- 
ductance dlo/dV vs. on-site energy eo and bias voltage V for 
A = hu v /3, A = hjj v , and (a) T = hw v and (b) F = 0.05 hw v . 
Features labeled "1" and "2" do not involve vibron emission 
or absorption, whereas those labeled "a"-"h" are due to this 
polaronic effect. All features are discussed in the text. 



III. RESULTS AND DISCUSSION 

A. Differential conductance of the dc Josephson 
current 

We now discuss the modification of the Josephson 
current due to the polaronic effect. In this study, we 
assume the two leads to be identical superconductors 
and take the amplitude of the order parameters to be 
A = |A^| = \Aji\. In Fig. 2] we first present density plots 
of the dc differential conductance dlo/dV as a function 
of the bias voltage V and the on-site energy eo, which 
in a break-junction setup could be controlled by a gate 
voltage. Different polaronic features are observed de- 
pending on the molecule-lead coupling T. For large T, 
dc transport is dominated by coherent tunneling across 
the molecule without requiring the energy of the electron 
to be aligned with the molecular level. For Fig. [2^a), 
we have chosen T = X = huj v . As noted above, in this 




eV = 2A 



FIG. 3: (Color online) Schematic representations of the two 
most important processes in the MJJ that involve neither vi- 
brons nor the molecular level. Left panel: single Andreev 
reflection at the right molecule-lead contact. Right panel: di- 
rect tunneling from the lower gap edge in the left lead to the 
upper gap edge in the right lead. The filled (empty) shapes 
represent the energy states of the superconducting leads below 
(above) the Fermi energy. Energy scales of the superconduct- 
ing gap A and the bias voltage eV are indicated. The arrows 
denote possible transitions of the electron or hole. 
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FIG. 4: (Color online) Schematic representation of various 
vibron-assisted tunneling and Andreev-reflection processes in 
the MJJ. The labels are the same as in Fig. [2] The DOS of 
both leads is modified due to the polaronic effect. The filled 
(empty) shapes represent the energy states of the supercon- 
ducting leads below (above) the Fermi energy. The molecular 
level eo (blue/gray) is broadened due to the coupling to the 
leads. Energy scales of the superconducting gap A, the vi- 
bron energy huj v , and the bias voltage eV are indicated. The 
arrows denote possible transitions of the electron or hole. 



regime the decoupling approximation is expected to un- 
derestimate the polaronic effect on transport but to show 
resonances at the correct voltagesPS The numerical re- 
sults thus provide a reasonable qualitative description to 
the vibron-assisted tunneling features. In Fig. [2]ja), we 
find several features at fixed bias voltage, insensitive to 
£ - The feature labeled by "1" at eV — hu> v /3 — A is due 
to single Andreev reflection without vibron emission or 
absorption. This process is illustrated by the left panel in 
Fig. [3] The feature labeled by "2" at eV = 2huj v /i = 2A 
is due to direct tunneling from the lower gap edge at one 
side to the upper gap edge at the other, illustrated by 
the right panel in Fig. [3] 

The features labeled by the letters "a"-"e" in Fig. 
[2]ja) involve vibrons. All these features and their repli- 
cas shifted by integer multiples of huj v can be explained 



by the onset of vibron-assisted Andreev reflections or 
coherent-tunneling processes sketched in Fig. [4j where 
the density of states (DOS) has been modified co mpared 
to the conventional picture of Andreev reflectio n 1 16 ! 22 ! 23 ! 
to account for the polaronic effect. For instance, the weak 
feature "a" located at uj v — ujj arises from the resonance 
of the vibron and Josephson frequencies and has been 
invoked by Marchenkov et aZP to interpret the observed 
over-the-gap structure. The process is sketched in Fig. 
|4][a), where an electron from a singular edge of the left 
lead undergoes one Andreev reflection at the right lead 
and arrives at a singular edge of the left lead. In this 
process, one vibron is emitted. 



On the other hand, for small T, where the decoupling 
approximation gives quantitatively reliable results, the 
tunneling processes are sensitive to the position of the 
molecular energy level. Pronounced features arise when 
singular edges of the DOS are aligned with the level £o- 
In Fig. [2jb), we show dlo/dV vs. V and £o for weak 
coupling, r = 0.05 huj v . As the non-resonant tunneling is 
strongly suppressed, the features seen in Figj2ja) become 
much weaker or are even invisible in Fig. |2Tb) . Instead 
we see vibron-induced features with peak bias voltages 
depending not only on lo v and A but also on £o- Two 
features, labeled by "f" and "g", satisfy deV/de = — 2. 
They are due to the alignment of the gap edges with 
the molecular level where resonant sequential tunneling 
through the molecule plays a dominant role. The cor- 
responding processes are depicted in Figs. Qf ) and (g). 
Note that the level positions are renormalized due to the 
electron-vibron coupling. One can see that feature "g" 
displays a sharp rise in dI /dV followed by a narrow re- 
gion of negative differential conductance. This is due to 
the onset of resonant tunneling at the gap edges with sin- 
gular DOSP Feature "f" instead shows a broader struc- 
ture without negative differential conductance, since the 
occupation of the final state is different, namely nearly 
empty for "f" and nearly full for "g". Moreover, we 
identify another pronounced feature "h" moving with 
£o, as well as one of its replicas. The underlying pic- 
ture is shown in Fig. [4jh) . An electron starting from 
the molecular level emits a vibron and is Andreev re- 
flected. One could say that the electron traverses the 
bias voltage V one and a half times. Accordingly, the 
feature has an unusual slope of deV/deo = —2/3. As- 
suming that the molecular level is shifted by a gate volt- 
age V g as £o = £q — eV g , we predict the distinctive slope 
dV/dVg = 2/3 of this feature in a bias- voltage/gate- 
voltage map. Interestingly, the feature is confined to the 
voltage range huj v < eV < huj v + 2A. This is due to the 
fact that for eV > tko v + 2A, electrons prefer direct tun- 
neling, while for eV < hu> v , the process is blocked due to 
the Pauli principle. 
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the polaronic effect: The anomalous (off-diagonal) self- 
energies in Eqs. ( 16 ) and ( 17 1 contain a factor (—1)'. This 
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FIG. 5: (Color online) Density plot of the ac differential con- 
ductance of (a) the dissipative Josephson current 7f and (b) 
the non-dissipative Josephson current if vs. the molecular 
energy level £o and the bias voltage V. The parameters are 
identical to those used for Fig. |2jb) . 



B. Differential conductance of the ac Josephson 
current 



We finally turn to the ac Josephson current. In Fig. 
[5] we plot the differential conductances of the dissipative 
and non-dissipative components as functions of V and 
£o for small T, where the decoupling approximation is 
valid. The features seen in Fig. |2jb) for the dc current 
are found again. However, their appearance is different: 
Feature "f" becomes blurred in the ac case and feature 
"h" is visible in a much broader voltage range. 

More interestingly, we observe an approximate an- 
tiperiodic behavior of both the dissipative and the non- 
dissipative components of the ac differential conductance 
as functions of the bias voltage. The ac current itself 
exhibits the same antiperiodicity (not shown). The an- 
tiperiod in eV is the vibron energy huj v . This means that 
the alternating current and the ac differential conduc- 
tance change their phase by 7r whenever eV is increased 
by HtUy. This antiperiodicity is a direct consequence of 



factor stems from the corresponding factor in the anoma- 
lous correlation function of the polaron-shift operators, 
Eq. (11). As noted above, the factor is due to Andreev 
reflection, since an outgoing electron and an Andreev- 
reffected hole couple to the vibron with opposite sign. 
For the weak molecule-lead coupling considered here, the 
current is dominated by processes involving a single vi- 
bron number I determined by V. According to Eqs. ( 16 ) 



and (|17j) , there is then an effective phase difference across 

The alternating current then 
As a result, the 



the MJJ of 2eVt/H + lw. 
contains a factor of sin(2eVt/h + lir). 
sign of the ac Josephson current depends on the even-odd 
parity of the vibron number I. The phase of the ac com- 
ponents can be measured with established techniques^ 
or employing the coupling to a charge qubit, as proposed 
recently.^ Ou r resu lts are related to the long-standing 
cos <p problem! 24 ! 25 ! It was found that the measured re- 
sult for the phase of the ac Josephson current does not 
agree with theoretical predictions. In the present work 
we have identified a mechanism by which this phase could 
even change periodically as a function of the bias voltage. 



IV. SUMMARY 

In summary, we have studied the transport properties 
of MJJs for which the electronic tunneling rates are mod- 
ified by polaron formation. Pronounced features due to 
its interplay with the superconductivity in the leads have 
been identified in the differential conductance of both the 
dc and ac Josephson currents. We have explained these 
features in terms of vibron-assisted Andreev reflection. 
The combination of sequential tunneling and Andreev re- 
flection leads to conductance peaks that show an unusual 
shift of their peak bias voltage with the molecular energy 
level or gate voltage, V ~ (2/3) V^. Furthermore, the 
opposite sign of the coupling of electrons and Andreev- 
reflected holes to vibrons induces periodic phase changes 
of the ac components of the Josephson current — their 
phase changes by 7r when the bias voltage eV is increased 
by one vibrational energy quantum huj v . We propose to 
search for this clear-cut polaronic effect by measuring the 
ac Josephson current through molecular junctions. 
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